Evaluation of the ^^Cr-^^Cr interaction via spin-flip scatterings 
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In order to evaluate go, the interaction strength of a pair of ^^Cr atoms with total spin 5 = 0, 
a specially designed s-wave scattering of the pair has been studied theoretically. Both the incident 
atom and the target atom trapped by a harmonic potential are polarized previously but in reverse 
directions. Due to spin-flip, the outgoing atom may have spin component jj, ranging from -3 to 3. 
The outgoing channels are classified by fi. The effect of go on the scattering amplitudes of each of 
these /^—channels has been predicted. 
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The chromium atoms ^^Cr are special due to having a 
larger spin F = 3 and a larger magnetic moment Gfis- 
Since the experimental realization of the chromium con- 
densation the interest in this field has been increas- 
ing. The ^^Cr condensate is a new kind of matter ag- 
gregation having the magnetic dipole-dipole interaction 
Vdd more than twenty times stronger than that of the 
alkalis family. A direct consequence of Vdd is the cou- 
pling of the spatial and spin degrees of freedom so that 
the conversion of spin angular momentum into orbital 
angular momentum can be realized. Thereby new phys- 
ical phenomena (say, creation of vortex) might appear. 
In addition to Vdd, the atom-atom interaction depends 
strongly on the total spin S of the pair, and is in gen- 
eral written as Vu = S{ri — J2s 9S^^ + ^dd, where 
gs is the strength and is the projection operator of 
the 5*— spin-channel (5 = 0, 2, 4 and 6). Up to now, 
02, fl4i and ge have been determined, but go has not yet 

0, Sj @|- However, many features of the condensate 
depend on this parameter (say, the phase-diagram of the 
ground state depends strongly on go 0, [^, 
dependence is also explicit in spin-evolution 
fore, the determination of go is important for a thorough 
and clear description of the condensates. This paper is 
one along this line. 

The effect of go can be exposed via two-body scatter- 
ings. We study the scattering of a very slowly incident 
^^Cr atom by another ^^Cr atom trapped previously in 
a potential U{r). The energy of collision is so low that 
the trapped atom is not able to be excited, and only the 
s-wave of the incoming one is involved. The scattering 
should be so designed that related observables are sensi- 
tive to go- The following points are proposed to meet the 
goal. 

(i) When Ms, the Z-component of the total spin S of 
the pair, is non-zero, go plays no role. Therefore, the 
two-body system must be designed to have Ms = 0. We 
therefore assume that the spin-state of the incident 
particle has fi — 3, while that of the target particle in the 
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FIG. 1; The s— wave cross section au of a particle scattered 
by a finite harmonic potential Up. ro is the width of the 
potential, and k = 0.1 is assumed. The units of au, ro, and k 
in this paper are A^, A, and A~^, respectively. A = y^h/mu! 
(say, 'du = 300 x 2tt, then A = 806nm for the ^^Cr atom). 



trap has fi = —3. This is different from the design of [1], 
where both particles have =-3. With their design, they 
have successfully determined the values of g2, gt, and ge. 
However, go can not be thereby determined. 

(ii) In general, for a scattering by a potential U{r), the 
cross-section depends on the parameters of the potential 
very sensitively. For an example, we study the case of a 
finite harmonic potential Up{r) — ^mLu'^{r'^ — Tq) if r < 
ro, or Up{r) = if r > tq. The parameter ro measures 
the width and depth of the potential. In what follows 
huj and A = ^h/mu! are used as units of energy and 
length, respectively. When a single ^^Cr atom with an 
incident kinetic energy is scattered by this potential, 
the s— wave cross-section ajj against the parameter ro 
is shown in Fig.l, where k = 0.1. One can see, if ro 
is given in a narrow domain around 2.512, au would be 
much larger. In the follows two kinds of potentials will be 
used, and their parameters will be selected in the narrow 
domains that lead to large ajj. 
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(iii) The dipole-dipole interaction is given by 

where is the operator of the spin of the z— th atom, 
and r — r2— ri. It can be proved that Vdd must cause an 
additional spatial relative rotation of the two particles in 
d-wave as shown in the appendix. However, if the to- 
tal orbital angular momentum L can be kept zero (say, 
due to a constraint in energy), the effect of Vdd would 
be strictly suppressed, and therefore can be neglected. 
The neglect of Vdd would greatly facilitate the calcula- 
tion shown below. This is different from the design of 
where a magnetic field is applied. It is reminded 
that, if higher partial waves are introduced by Vdd, they 
are accompanied by large-S" components due to the con- 
servation of the total angular momentum. Among these 
components those with a large- Ms will lead to a decrease 
of energy due to the Zeeman effect. Therefore, the in- 
crease of rotation energy can be partially canceled by 
the additional Zeeman energy. Consequently, the mag- 
netic field would help Vdd to cause a stronger mixing to 
include higher partial waves, and the effect of Vdd thereby 
becomes remarkable. This is the case of 

With the above consideration, for numerical calcula- 
tions, two potentials are introduced. One is the finite har- 
monic well Up(r) with the parameter tq mentioned above, 
the other is the Gaussian potential Ua{r) = — 4e~'^''/''^^ 
with the parameter a. The propagation number of the 
incoming atom k < 0.2 so that the related kinetic en- 
ergy is much smaller than the energy gap of the 
ground state, and the maximal effective classical angular 
momentum of the incident particle hkro is smaller than 
h. With these choices the target atom can not be excited, 
and only the s-waves of both atoms are involved, and the 
total orbital angular momentum L remains to be zero as 
required. 

When Vdd is dropped, S and Ms are good quantum 
numbers. Let the total wave function ^E* be a composition 
of different S"— component, 

* = VC5*5 = VCs — $5('^l,r2)(12)sMs (2) 
s s ''""^ 

where S is from to 6 and Ms = as designed. Each s 
should satisfy the Schrodinger equation, while the coef- 
ficients Cs are determined by how the incoming channel 
would be. Since only s-waves are involved, the angular 
degrees of freedom are irrelevant. Thus the Schrodinger 
equation for <^s reads 

-5s5(ri-r2)] $s(ri,r2) = (3) 

Where the total energy E = + fc^/2, and Eg is the 
ground state energy (negative) of the target particle. In- 
cidentally, Eg and the associated wave function (f>g are 



easy to be obtained. In order to facilitate numerical cal- 
culation, the S — function is replaced by e~^^ /(tt//?)'^/^, 
where /? is chosen to be a very large number to assure a 
very short range. We choose f3 = 1/(0.06)^, and we found 
that the effect of the change of /3 in a reasonable domain 
is slight. When S is even, $s must be symmetric with 
respect to the interchange of ri and r2 , whereas when S 
is odd, it must be anti-symmetric. 

For fc = 0.1 and U — Up, we know from Fig.l that 
the optimal tq leading to the largest ajj is 2.512. Ac- 
cordingly, three cases with ro =2.5, 2.512, and 2.53 are 
chosen to see the effect of the variation of rg around its 
optimal value. We define a domain (0,5) for both ri and 
r2. Then, Eq.(3) is solved numerically in the domain. 
The crucial point is the embedment of boundary condi- 
tions. Obviously, we have ^s{Q,i'2) = $5(^1,0) — 0. 
Furthermore, it is reminded that, due to the constraint 
in energy, if one atom is far away, the other one would 
remain in the trap in the ground state (the outgoing 
of both atoms is not possible). Therefore, when r2 is 
sufficiently large, U{r2) = 0, V12 — 0, and $s(ri,r2) 
would tend to (j)g{ri){bs cos{kr2) + 05 sin(fcr'2)), where 
the second factor is the common form of s-waves. Both 
6s and as are real numbers, and the unitarity is thereby 
assured. When r2 = 5, the second factor is just a con- 
stant denoted by hs- Now we arrive at the second pair of 
boundary conditions, namely, <I>s(ri,5) = 4>g{ri)hs and 
*i's(5,?'2) = {—lY (f)g[r2)hs- Where, hs can be arbitrary 
given because the eventual results do not depend on it as 
shown below. The factor (—1)"^ is needed to assure the 
correct permutation symmetry. 

Now, the values of $s('"i,''2) at all the borders are 
specified. Thus the problem is to solve an elliptic par- 
tial differential equation under the Dirichlet condition. 
Numerical solutions can therefore be obtained by using 
standard programs of difference equations. For S = 2, 4, 
and 6, gs have been known. Using Tiuj and A as units, 
92 = -0.134 X 10-3^, 54 = 1.11 X 10-3^, and 
36 = 2.14 X lO^^V^. For S ^ 0, gs will be given at 
a number of values from —2(76 to 2g(,. When S is odd, gs 
is simply zero. Examples of ^s are given in Fig. 2. One 
can see in Fig. 2a that the two atoms can be very close 
to each other in the trap if S is even, but can not if S is 
odd as shown in 2b. Obviously, <I's(''ij^2) is symmetric 
(anti-symmetric) in 2a (2b) as expected. 

Once $s has been solved numerically, we can extract 
hs and as from its asymptotic form. It is reminded that 
the incoming particle has = 3, and the target particle 
has /i = — 3. To meet this requirement, we choose Cs = 
^3 3 3 -s/iO'S — ibs)- With this choice, in the asymptotic 
region (say, ri is large while r2 is small), the total wave 
function appears as 

* ^ — 0g(r2)[sin(fcri)x3(l)X-3(2) 

+ ^^exp(zfcri)xp(l)X-p(2)] (4) 
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FIG. 2: Contour plots of $s{ri,r2) with S = 6 (a) and 5" = 
odd (b). The potential is Up{r) with ro = 2.5, and k = 0.1. 
uj = 300 X 2ti is assumed (the same in the following figures). 

where 

where both even and odd S are included in the sum- 
mation, is the s-wave scattering amplitude of the 
//—channel (i.e., the outgoing particle has spin compo- 
nent /i, while the inner one has — /i ). The asymptotic 
form with a large r2 is similar. From (4) we know that, 
in addition to the scattering by the potential, spin-flips 
might occur due to the spin-dependent interaction V12. 
Incidentally, since Cs contains the factor (05 — ibs)~^ 
and both as and 65 are proportional to hs, depends 
not at all on the set of constants hs chosen before. 

We are interested in how the s-wave cross section of 
a /i— channel ct^ = 47r|/^p/A:^ is affected by go- It was 
found that (73 (elastic channel without spin-flip) is much 
larger than other tr^. When the cross section of pure 
potential-scattering au is considered as a unit, the ratio 
a-ijou against go is plotted in Fig. 3. Four points are 
noted: 

(i) The magnitude of (T3 is in general close to aij. In 
particular, if go is close to gi, ""a would be very close to 

CFU- 

(ii) (T3 depends sensitively on the potential parameter, 
because uu does. It implies that an appropriate choice 
of parameters is crucial. 

(iii) When the potential is narrower than the one lead- 
ing to the optimal au^ if go is smaller than g2, a more 
negative go would lead to a larger 0-3, while a more pos- 
itive go lead to a smaller <jj, (refer to the dash curves 
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FIG. 3: (color online), a-i/au against go/ga with k = 0.1. 
Finite harmonic potential Up and Gaussian potential Ua are 
used in a and b, respectively. In 3a, the parameter ro is 
given at 2.5, 2.512, and 2.53 marked by the curves. The 
corresponding au are 814, 1257(the optimal value), and 779, 
respectively. In 3b, the parameter a is given at 1.48, 1.50, 
and 1.52. The corresponding au are 978, 1257 (the optimal 
value), and 1003, respectively. 



of Fig. 3). When the potential is broader, the effect of 
go is in reverse (refer to the dash-dot-dot curves). How- 
ever, when the parameter of the potential is optimized, 
a larger |go — g2| would lead to a smaller (73 (refer to the 
bold curves). 

(iv) The above qualitative feature does not depend on 
the details of the potential. This is explicit by comparing 
3a and 3b, where distinct potentials are used. Thus, for 
selecting an appropriate potential, the crucial point is the 
understanding of the optimal parameter. 

Incidentally, in all the cases of Fig. 3 the ground state is 
deeply bound. For an example, when Up{r) has rg = 2.5, 
the ground state energy is -1.6294, while the first excited 
state is about 1 higher. Thus the gap is greatly larger 
than the bombarding energy fc^/2 < 0.02. 

In Fig. 4 17^/(73 against go is plotted. Three points are 
noted: 

(i) The cross-sections of spin-flip channels are re- 
markably smaller than that of the non-spin-flip channel. 

(ii) It is reminded that the three with odd S are 
free from atom-atom interaction, therefore they have 
the same hs/{as - ihs). Since E's C'^s, 3,-3^3;° 3,-^ = 
(<5/i,3 — <^M,-3)/2, where the summation covers only the 
three odd 5, it is straight forward to prove = if 
/i ^ ±3. Hence, the cases with /i = — 2 and -1 can not 
be seen in the figure. 
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FIG. 4: (color online), (J^/a-i against go/ ga. The potential Up 
with ro = 2.512 (the optimal value) is used, and k — 0.1. n is 
marked by the curves, a-i overlaps ai, and a-2 overlaps (T2. 
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FIG. 5: (color online), The same as Fig. 3a but with k — 0.2. 

(ill) Besides the non-spin-flip channel, the most impor- 
tant spin-flip channels would be the fi = ±2 channels if 
go is negative. However, if go is positive, /i = chan- 
nel would be a little more important, while a±2 is the 
smallest. In particular, when go 56/2, a±2 is close to 
zero. 

The features of Fig. 4 are quite popular. When dif- 
ferent potentials with different parameters are used, the 
qualitative features remain unchanged. 

In previous examples fc = 0.1 is taken. When k varies. 



the above qualitative features remain. This is shown in 
Fig. 5 to be compared with Fig. 3a. Note that, when k 
varies, the optimal parameters of the potentials would 
vary accordingly. 

In summary, a specific s-wave scattering of a pair of 
^^Cr atoms has been designed and studied theoretically. 
Different from the design of , both atoms are polarized 
but in reverse directions previously and no magnetic field 
is applied. The effect of go on the scattering amplitudes 
without and with spin-flip has been predicted. Whether 
go is smaller or larger than 52 is found to be crucial to the 
qualitative features. In order to determine the unknown 
go, associated experiments are desirable. 
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APPENDIX A: MATRIX ELEMENTS OF THE 
DIPOLE-DIPOLE INTERACTION 

In this appendix the matrix elements of the dipole- 
dipole interaction between the total spin-states (12)5^ 
of a pair of spin-3 atoms are given. 

Making use of the spherical components of the spin- 
operator Fi and F2, we have 

((12)s-M' \Vdd\il2)sM) = -^252V4^Ci%° ^^0^23^ 
r 1 1 2 1 J 

< 3 3 > C'2.M-M'. S'.Af'^2,M-M'(f) (Al) 

[ 3 3 S" J 

where the Clebsch-Gordan and 9-j symbols fld\ are in- 
troduced, and S + S' must be even. From this formula 
it is clear that Vdd must cause a spatial relative rotation 
in d-wave together with a change of S (the choices of S' 
is |S'-2|, S, and 5-^2). 
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